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Abstract -We consider a solution of the effective four-dimensional brane-world equations, 
obtained from the General Relativistic Schwarzschild metric via the principle of Minimal 
Geometric Deformation, and investigate the corresponding signatures stemming from the 
possible existence of a warped extra-dimension. In particular, we derive bounds on an extra- 
dimensional parameter, closely related with the fundamental gravitational length, from the 
experimental results of the classical tests of General Relativity in the Solar system. 


Introduction. — Brane-world (BW) models [T| 
represent a well-known branch of contemporary high- 
energy physics, inspired and supported by string the¬ 
ory. These models are indeed a straightforward 5D phe¬ 
nomenological realisation of the Horava-Witten super¬ 
gravity solutions [2], when the hidden brane is moved 
to infinity along one extra-dimension, and the moduli 
effects from the remaining compact extra-dimensions 
may be neglected 0- The brane self-gravity, encoded 
in the brane tension er, is one of the fundamental 
parameters appearing in all BW models, with cr -1 / 2 
playing the role of the (5D) fundamental gravitational 
length scale @. In this work, we shall explicitly study 
the observational effects determined by the parame¬ 
ter /3 ~ (cr -1 / 2 /!?) 2 , which describes a candidate for 
the modified 4D geometry surrounding a star of ra¬ 
dius i?in the BW. This particular geometry will be 
obtained as an exact minimal geometric deformation 

1 We shall mostly use units with the four-dimensional Newton 

constant G = c = 1, unless otherwise specified. 


(MGD) [J] of the Schwarzschild solution to the field 
equations in General Relativity (GR). The MGD ap¬ 
proach ensures, by construction, that this BW solution 
smoothly reduces to the GR Schwarzschild metric in 
the limit er -1 — > 0, thus allowing us to analyse varia¬ 
tions from GR predictions for small values of the de¬ 
forming parameter /?. This parameter controls the cor¬ 
rections and is related to the brane tension, the radius 
of the star, and will also be shown to depend on the 
compactness of the star. 

Minimal Geometric Deformation. — The effec¬ 
tive Einstein equations on the brane take the form j5] 

Gfj,v — Ijj u A.g^ , ( 1 ) 

where T^ u = T^ u +1 S^ + gl- £^ denotes the effective 
energy-momentum tensor, with the stress tensor of 
brane matter, and and S ^ the non-local and high- 

energy Kaluza-Klein corrections. If BW matter is a 
perfect fluid with 4-velocity id*, and h^ = g^ 
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are the components of the metric tensor orthogonal to 
the fluid lines, then 


I(r, n>) 




2i/ 

x 

2 

x 



( 6 ) 


Spv = 


IA [ u^ Ui/ T ^ 


+ 'P f iv J rQ( l i Mj/) 


, ( 2 ) 


where U denotes the bulk Weyl scalar, V^ is the 
anisotropic stress and Q M the energy flux. 

Solving the effective 4D Einstein equations in the 
BW is a hard task and, already in the simple case of a 
spherically symmetric metric, 


ds 2 = e v dt 2 — e x dr 2 — r 2 dQ 2 , (3) 


only a few “vacuum” solutions are known analyti¬ 
cally ©no]. Moreover, for static matter distributions 
Qfj, = 0 and = V (r M r v + ± h^) 0. For stellar 
systems, the quest for BW solutions becomes even more 
intricate, mainly due to the presence of extra terms, 
non-linear in the matter fields, which emerge from high- 
energy corrections [51I51HB], Nonetheless, two approxi¬ 
mate analytical solutions have been found in the MGD 
approach. It is worth to emphasize that these met¬ 
rics are exact solutions of the effective equations ©, 
although they are not complete solutions of the full 
5D equations mm- This approach also yields physi¬ 
cally acceptable interior solutions for stars [T9] , relates 
the exterior tidal charge found in Ref. [TO] to the ADM 
mass, and let us study (micro) black holes [20I2T1 . eluci¬ 
dates the role of exterior Weyl stresses from bulk gravi¬ 
tons on compact stellar distributions [22] and shows 
the existence of BW stars with Schwarzschild exterior 
without energy leaking into the bulk [23]. Moreover, 
both the associated 5D solutions and black strings were 
obtained in various contexts jT21ll5ll24ll25] . with mod¬ 
els for the quasar luminosity variation induced by BW 
effects MMH- 

Let us start by revisiting the MGD approach, which 
is built on the requirement that GR must be recovered 
in the low energy limit cr _1 —> 0. In particular, by 
solving the effective 4D equations ©, the radial com¬ 
ponent of the metric is deformed by bulk effects and 
can be written as [21. 


e A = p + / , 


where 



(4) 

(5) 


and 



1 - 


2 m[r ) 
r 


for r > R , 
for r < R , 


where m denotes the standard GR interior mass func¬ 
tion. The constant M depends in general on the brane 
tension a and must take the value of the GR mass 
Mo = m{R) in the absence of BW effects, namely 
M 0 = M| ct -i = 0 . The function H in © is given by 


H(P,P, v) 


24 7r p — [j! I — H— 


v’ 1 

2 r 

n 2z/ 


-\-p I v -|——-|-1—o 


1 

r 2 


,(7) 


and encodes anisotropic effects due to the bulk grav¬ 
ity on the pressure p , matter density p and the metric 
function v. Finally, the parameter jd in © depends 
on the brane tension cr, the radius R and the mass 
M of the self-gravitating system, and must be zero in 
the GR limit. In the interior, r < R, the condition 
/3 = /3(cr, R, M) = 0 must hold in order to avoid singu¬ 
lar solutions at r = 0 [since the integral in © would 
diverge for ro —> 0]. However, for a vacuum solution, or 
more properly, in the region r > R where there is only a 
Weyl fluid surrounding the spherically symmetric star, 
the parameter (3 can differ from zero. 

The crucial point is that, any given perfect fluid 
solution in GR yields H(p,p,v) = 0, which provides 
the foundation for the MGD approach. In fact, every 
perfect fluid solution in GR can be used to produce 
a minimal deformation on the radial metric compo¬ 
nent in the sense that all the deforming terms in 
Eq. © are removed, except for (a) those produced 
by the density p and pressure p, which are always 
present in a realistic stellar interior (where (3 = 0 
for r < R), and (b) the one proportional to the pa¬ 
rameter /3 in a vacuum exterior (with p = p = 0 
for r > R). It is worth to emphasise that the con¬ 
dition H = 0 holds for any BW solution obtained 
by the MGD approach, and corresponds to a mini¬ 
mal deformation in the sense explained above. More¬ 
over, H may not be negative when a perfect fluid is 
used as the gravitational source on the brane, since [4] 
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H(p, p , u) = 24 7T p— (2 G 2 2 + G \) |o—i = o and the com¬ 
ponents of the Einstein tensor G 1 1 = G 2 2 = 8 7rp for a 
spherically symmetric perfect fluid, so that the condi¬ 
tion H = 0 always holds. 

In order to obtain a deformed exterior geometry, 
we then start by inserting the spherically symmetric 
Schwarzscliild metric 


e vs = e~ Xs = 1 - 


2 M 


( 8 ) 


in the expression a for r > R, where p = p = 0. 
Since © is a GR solution, H(r > R) = 0 and the 
correction in Eq. © will thus be minimal, 


/+(r) = f{r)\ p -p- H - 0 — P e 


-I 


(9) 


The outer radial metric component © will read 

j+ 2 M j 

e~ x = 1-+ (3(a, R, M) e~ J , (10) 

r 

which clearly represents a BW solution different from 
the GR Schwarzschild metric, with /3 equal to the 
extra-dimensional correction to the GR vacuum evalu¬ 
ated at the star surface, that is /3 = / + (r = R), and 
containing a “Weyl fluid” for r > R [22] . 

We next consider the general matching conditions 
between a general interior MGD metric (for r < R), 

rlr 2 

ds 2 = e v {r) dt 2 -— n - r 2 dfl 2 , (11) 

1 -*2fel+/-(r) 


where / is also given by Eq. © with H = 0, and 
the above exterior metric (for r > R), which can be 
written like m by replacing — with +. Continuity of 
the metric at the star surface E of radius r = R yields 


'Ri 


2 M 

~R 


2 M 0 
R 


+ {fn ~ fR) , ( 12 ) 


where F R = F(r R±) for any function F. Continu¬ 
ity of the second fundamental form on E likewise pro¬ 
vides the expression [G p!/ r"] s = 0, where denotes a 
unit radial vector and [/]e = f(r —> R + )~ f(r —► R~). 
On using the effective 4D equations ©, this condition 
becomes 


2 U 


EL 

2 


pp 


V 


= 0 . 


(13) 


J E 


Since the star is assumed to be only surrounded by a 
Weyl fluid described by U + and V + (and p = p = 0) 
for r > R, this matching condition takes the final form 

ap R +4F R +2U R + ^ + PrPr = 2(2V+ +U+) (14) 


with p R = p R and p R = p R . The limit <r _1 —> 0 in 
Eq. (1141) leads to the well-known GR matching condi¬ 
tion pr = 0 at the star surface. Eqs. and (ED are 
the necessary and sufficient conditions for the matching 
of the interior MGD metric to a spherically symmetric 
“vacuum” filled by a BW Weyl fluid [28]. 

BW star. — BW effects on spherically symmetric 
stellar systems have already been extensively studied 
(see, e.g. Refs. ©J] for some recent results). Let us 
now investigate in details the MGD function f + (r) pro¬ 
duced by the Schwarzschild solution ©. By inserting 
it into Eq. ©, we obtain 


f+(r)=p(a,R,M)- 

r 


1 _ 2M 


1 3 M ’ 

1 2 r 


(15) 


where b is a length given by b = R( 1 — |^)/(1 — 2M) 
and the deformed exterior metric components read 




2 M 



/?(er, R, M) b 

1 + 1 _ 3M „ 

L 2 r 


(16a) 

(16b) 


matching the vacuum solution found in Ref. [25] when 
j3b = K/a, with K > 0. The corresponding Weyl fluid 
is described by [221 

v+ p b ( 1-tJf) U+ —/3bM 

a 9r s(l-*M) 2 ’ ^ 12r .4 (l _m) 2 


We can now obtain the parameter /? = j3(a,R,M), 
depending on the interior structure, by employing the 
deformed Schwarzschild metric (I Rial) and (Il6bl) in the 
matching conditions (fI21) and ED- Eq. 1121) just be¬ 
comes e v R = 1 — , whereas Eq. (fTH) yields 

VR+ij k^' R + ll) = -%' (18) 

with v' R = {y~y \ r —R- These are the necessary and 
sufficient conditions for matching the two minimally 
deformed metrics given by Eqs. (fill) . (I16al) and (I16bl) . 
If M in Eq. © were the GR mass Mo, one would have 
f R = f R [see Eq. (1 1 2[) ] . which is an unphysical condi¬ 
tion, according to Eq. ([TSf) . In fact, the interior defor¬ 
mation / = f—(r) is positive, but the matching condi¬ 
tion (fTHl) shows that the exterior deformation must be 
negative at the star surface, f R < 0, or else a negative 


P-3 





















R. Casadio et al. 


pressure pr < 0 would appear. Hence, according to 
Eq. (USD, the deformation / + (r > R) is negative for 

/3 < o E3j. 

The exterior geometry given by Eqs. (11 6 all and (116bl) 
may seem to have two horizons, namely rh = 2 M and 
r 2 = 3 M/2 — fib. However, since /3 must be propor¬ 
tional to cr^ 1 in order to recover GR, the condition 
r 2 < Th must hold, and the outer horizon radius is 
given by r/j = 2 M. The specific value /3 = —M/ 2 
would produce a single horizon = r 2 = 2 M, but the 
limit er _1 — > 0 does not reproduce the Schwarzschild so¬ 
lution, as seen from the condition Mq = M | CT -i =0 . On 
the other hand, f R < 0 implies that the deformed hori¬ 
zon radius = 2 M is smaller than the Schwarzschild 
radius r# = 2 M 0 , as it can be clearly realised from 
Eq. (HU). This general result shows that 5D effects 
weaken the strength of the gravitational field produced 
by the self-gravitating stellar system. 

Finally, when (1151) is considered in the matching con¬ 
dition (USD , we obtain 


0 = it = -R 2 


PR + 


[r 

R 


(19) 


showing that /3 is always negative and (interior) model- 
dependent through v' R . In particular, we can find /3 by 
considering the exact interior BW solution of Ref. El. 
that is 


r 


32 C \ 240+589C'r 2 -25C' 2 r 4 -41C' 3 r 6 -3C' 4 r 8 
49 ct [ 3(1 + Cr 2 ) 4 (l + 3CY 2 ) 


80 arctan(-\/CV) 
(1 + CY 2 ) 2 (1 + 3 Cr 2 )VCr 


( 20 ) 


where C denotes a constant (with the same dimen¬ 
sions of er) given by C R 2 = = a, and v' = 

8Cr (1 + CY 2 ) -1 . Using the explicit form of f~(R) and 
Pr = 0 in Eq. (fl9D yield 


/?0, R) = f+ = 


Co 

R? a ' 


( 21 ) 


where Co — 1-35 is a (dimensionless) constant. The ex¬ 
terior deformation is finally obtained by using Eq. m 
in Eq. (USD, leading to 


f + = 


C 0 b / r — 2 M 0 \ 
R 2 a r \ r — 3 M 0 ) 


+ C(a~ 2 ) , 


( 22 ) 


monotonically increasing function of the star compact¬ 
ness Mq/R. Since extra-dimensional effects are the 
strongest at the surface r = R and become more im¬ 
portant for smaller stellar distributions, the more com¬ 
pact the star the larger /?, and thus the MGD of the 
GR solution. 

Solar System Classical Tests. — Classical tests 
in the Solar system can probe BW signatures. The per¬ 
ihelion precession of Mercury, the deflection of light by 
the Sun and the radar echo delay observations are well- 
known tests for the Schwarzschild solution of GR and, 
in BW models, for the DMPR and the Casadio-Fabbri- 
Mazzacurati metrics as well. BW effects in spherically 
symmetric space-times were comprehensively studied, 
e.g., in Ref. [30]. In our case, Solar system tests will be 
employed to bound the MGD parameter /3 in Eq. m- 

Perihelion Precession. A test particle in a spheri¬ 
cally symmetric metric (0 has two constants of motion, 
E and L, respectively yielding energy and angular mo¬ 
mentum conservation. By the usual change of variable 
r = 1/u and defining 

g(u) = 1 - e" A , (23) 


the relevant equation of motion reads [303 


cPu _ 1 d / E 2 e~ x ~ u 
d(f) 2 + U 2 du\ c 2 L 2 


— A 


-jY+g{u)w 


= k(u) . (24) 


1/2 

By denoting y(u) = (1 — ( dk/du ) | Uo ) ' , a circular or¬ 
bit m = Mo is determined by the root of the equation 
uq = k(uo), and a deviation with respect to it is pro¬ 
vided by S = 6q cos + a), with do and a con¬ 

stants [3U]. The variation of the orbital angle with re¬ 
spect to successive perihelia is given by cj> = , 

where the perihelion advance is i ~ ^ u - u ’ ^ or 
small values of ( dk/du) u=Uo . For a complete rotation, 
the perihelion advance is S(f> ~ 2ttl. 

We now consider the perihelion precession of a planet 
in the MGD geometry described by Eqs. (llfial) and 
(llfibl) . Since L is related to the orbit parameters by 
L = 27ra 2 v / l — e 2 /cT [30], where T denotes the period 
of motion, the perihelion advance thus yields 


8<t> = 5(i>GR ~ f(/3) , (25) 


where bo = b(Mo) is given by the length b at M = where S(f>GR = 67 iGM/c 2 a (l — e 2 ) is the well-known 
Mo- The deformation / + (r > R) is therefore a Schwarzschild precession formula and f(/3) ~ 673.94/3, 
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where we employed c = 2.998 x 10 s m/s, Mq = 1.989 x 
10 30 kg, a = 57.91 x 10 9 m, R 0 = 6.955 x 10 8 m, 
e = 0.205615, and G = 6.67 x 10 -11 m 3 kg -1 s -2 . 

The observed difference Scj> — ScpoR = 0.13 ± 
0.21 arcsec/century [Jl] can thus be ascribed to BW ef¬ 
fects, according to Eq. (|25j). Observational data Eqeb 
yield the bound f(f3) < (1.89±2.33) x 10 -8 , which con¬ 
strains 

/3 < (2.80 ±3.45) x 10" 11 . (26) 

Light Deflection. A similar procedure describes 
photons on null geodesics, with the equation of motion 
that can be written as 

(^) + u2 = ^j 2 e ~ V ~ X + 9(u)u 2 = p{u) , (27) 


which therefore implies + u = ^ dp }™' > ■ In the lowest 
approximation, the solution is u = c ^ ^ , where R 0 is 
the distance of closest approach to the mass M. It can 
be iteratively employed in the above equation, yielding 


COS (f> ' 

Rq 


^ jdu. The total deflection 


cPu/dcj) 2 + u = p ^ 
angle of the light ray is given by S = 2e 
For the geometry (1 16 all and (Il6bl) . Eq. 
g(u) = (2 GM/c 2 ) u, resulting in 


leads to 


p(u) = 


flbo 


(2-3 OMu^ 


GMu ( 9GM 


0 


-u-11 +3 


- 2a z ) + 


3 GMu 2 


(28) 


since r = V x 2 + R 2 - The above integrand takes the 
form 

»p(A0- 1 = (i-^r/snfev+i)’' 

„ 2 GM flip 4pe 0 GM 

c 2 r 3GM — 2c 2 r c 2 r(3GM—2c 2 r) ’ 

where we used a first order approximation based on 
Eqs. (I16al) and (I16bl) . Eq. (1311) . using the approxima¬ 
tions R 2 /i 2 1 (* = 1,2), and considering terms up 

to order (GM/c 2 R) 2 , hence reads 


ST - ST gr + 


c 3 R 


In 


R 2 


5t tGM 


(33) 


which reproduces the Schwarzschild radar delay 
<5 Tgr = 2( / 3 M hi when (3 = 0, and the second 
term imposes a constraint on BW models. Recent 
measurements of the frequency shift of radio photons 
both to and from the Cassini spacecraft, as they passed 
near the Sun, have refined the observational constraints 
on the radio echo delay. For the time delay of the 
signals emitted on Earth towards the Sun, one ob¬ 
tains Af radar = A t (1 T A radar ), with A radar — 

(1.1 ± 1.2) x 10 -5 [32] - In the BW geometry (U6al) 
and (Il6bl) . measurements of the frequency shift of ra¬ 
dio photons [30l[32] finally yield the physical bound 


57rG 2 M 2 A radar 

2£ 0 R\n(^) 


(3.96 ± 4.30) x 10 -5 . (34) 


The total deflection of light is finally given by 


5(/) 


3GM 


+ /3 b 0 


E 2 Rq 

c 2 L 2 


1871C 2 Rq \ 

GM ) ’ 


(29) 


in the limit ( , ( 1 Ro j « 1, f C 1 and -p— 1 <C 1, 
which implies the bound 


/3 < (1.07 ±4.28) x 1CT 10 . (30) 


Radar Echo Delay. Another classical test of GR 
measures the time for radar signals to travel to, for 
instance, a planet [3T|. The time for light to travel 
between two planets, respectively at a distance l\ and 
^2 from the Sun, is well-known to be To = / i dx/c. 
On the other hand, if light travels in the vicinity of the 
Sun, the time lapse (ST = T — T 0 is given by [30] 



(31) 


This provides a bound on the MGD parameter /3, which 
is the weakest one among those in our analysis. 

Concluding Remarks. — BW models can be con¬ 
fronted with astronomical and astrophysical observa¬ 
tions at the Solar system scale. In this paper we have 
in particular considered the BW exterior solution (U6al) 
and (I16bl) obtained by means of the MGD procedure, 
and compared its predictions with standard GR results. 
This exterior geometric contains a parameter f3 and we 
were able to constrain it from the presently available 
observational data in the Solar system. We found the 
strongest constraint is given by measurements of the 
perihelion precession, namely Eq. (1261) . 

Let us recall that limits for the brane tension in the 
DMPR and Casadio-Fabbri-Mazzacurati BW solutions 
have already been determined via the classical tests 
of GR [3D\ Since bounds on the parameter fl imply 
lower bounds for the brane tension from Eq. m , we 
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can conclude that the constraint (1551) complies with 
the ones provided by such solutions of the effective 4D 
Einstein equations ©■ In fact, the brane tension in 
the MGD framework is bounded according to 


9 M q c 2 

(i 2GA M 2 

V C 2 Rq ) 

kRq/3 

(, 3GM e \ 

V 2 


(35) 


which implies that a > 5.19 x 10 6 MeV 4 , when the 
bound (1261) is taken into account (we omit errors here 
since we are solely interested in orders of magnitude). 
This bound is still much stronger than the cosmological 
nucleosynthesis constraint, however much weaker than 
the lower bound obtained from measurements of the 
Newton law at short scales. We can therefore conclude 
that the MGD geometry (U6aD and (116bl) is acceptable 
within the present measurements of BW high-energy 
corrections. 

After this work had been completed, we developed 
an extension of the MGD, which produces a rich but 
complex set of new exterior solutions [35], whose com¬ 
plete analysis is highly non-trivial. However, the solu¬ 
tion used here represents the simplest non-trivial exten¬ 
sion of the Schwarzschild solution within the extended 
MGD approach. 
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